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The energy transfer processes and coherent phenomena in the fucoxanthin-chlorophyll protein complex, which 
is responsible for the light harvesting function in marine algae diatoms, were investigated at 77 K by using two- 
dimensional electronic spectroscopy. Experiments performed on the femtosecond and picosecond timescales 
led to separation of spectral dynamics, witnessing evolutions of coherence and population states of the system 
in the spectral region of Q y transitions of chlorophylls a and c. Analysis of the coherence dynamics allowed 
us to identify chlorophyll (Chi) a and fucoxanthin intramolecular vibrations dominating over the first few 
picoseconds. Closer inspection of the spectral region of the Q y transition of Chi c revealed previously not 
identified mutually non-interacting chlorophyll c states participating in femtosecond or picosecond energy 
transfer to the Chi a molecules. Consideration of separated coherent and incoherent dynamics allowed us to 
hypothesize the vibrations-assisted coherent energy transfer between Chi c and Chi a and the overall spatial 
arrangement of chlorophyll molecules. 
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I. INTRODUCTION 

Diatoms are unicellular chromophyte algae inhabit¬ 
ing marine environment. They are the major players in 
photosynthesis, accounting for nearly a quarter of the 
global primary production 1 3 . Light-harvesting in di¬ 
atoms is performed by their intrinsic membrane proteins, 
such as the fucoxanthin-chlorophyll protein (FCP) com¬ 
plex. FCP shares some structural and sequence homol¬ 
ogy with the light-harvesting complex II (LHCII) from 
the higher plants, but contains more carotenoids per 
chlorophyll. Carotenoids fucoxanthins (Fx) strongly ab¬ 
sorb blue-green photons with the longest propagation 
distances in water. Such adaptation allows diatoms to 
efficiently harvest sunlight in underwater conditions. 

Information on the molecular structure of FCP is very 
limited; there is no crystallographic data of the com¬ 
plex. Several preliminary models of pigment organiza¬ 
tion in FCP, based on the assumption that the structure 
of FCP is similar to the LHCII, have been proposed 3 ®. 
These models, however, are currently limited for sev¬ 
eral reasons. First, only some pigment binding sites are 
conserved between FCP and LHCII. Second, the pig¬ 
ment composition and, thus, detailed functional roles 
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in LHCII and FCP are markedly different, as the latter 
contain nearly as many fucoxanthin molecules as chloro¬ 
phylls (Chi)—most probably 8 Fx, 8 Chi a and 2 Chi c 
molecule^®— while LHCII is composed of 8 Chi a, 6 
Chi b and 4 xanthophyll molecules^*. 

Chi c molecules enhance the absorption of blue photons 
in FCP, as welffl. Ultrafast molecular excitation trans¬ 
fer from Chi c to Chi a molecules was shown to occur at 
room temperature 12 . Nonetheless, the whole picture of 
the energy funnels in the FCP complex and the physical 
phenomena ensuring ultrafast energy transfer still needs 
to be detailed. Unfortunately, lack of structural informa¬ 
tion about FCP prevents detailed numerical studies. All 
the information about the inner workings of this complex 
comes solely from the experimental data and its analysis, 
which is often done in the light of the proposed models^®. 
Therefore, it is important to obtain as much experimen¬ 
tal data as possible to constrain and refine the current 
models. 

The two dimensional (2D) electronic spectroscopy (ES) 
is a remarkable tool for the analysis of coherent and in¬ 
coherent phenomena appearing in molecular excitation 
dynamics. The first and the most recognized results of 
2D ES demonstrating complex behavior of the excitation 
evolution were delivered for the photosynthetic Fenna- 
Matthews-Olson (FMO) complex from purple bacteria 13 . 
Since then, 2D ES has been successfully applied for 
studies of a number of different systems; 2D ES es- 
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sentially helped to establish or update existing kinetic 
schemes of light-harvesting antennae and reaction cen¬ 
ters of plantd^ 17 and bacteriaP^S The 2D ES was also 
performed on the whole photosystem I complex demon¬ 
strating signatures of fast energy transfer 21 . In addi¬ 
tion, it was utilized in studies of other molecular aggre¬ 
gates, for instance, polymer j^. cylindrical (bi-tubular) 
J-aggregated^, quantum dotd^ etc. 

Most of these findings were obtained from experi¬ 
ments and theoretical simulations performed using (or 
assuming) the knowledge of the structural organization 
of the pigment molecules (often with resolution of a few 
angstrms). Thus, these studies were valuable not only for 
widening a general understanding of the physical mech¬ 
anisms within molecular systems, but also were conve¬ 
nient for testing and interpreting the outcomes of the 
new method of 2D ES. As the method itself has already 
been well approved and appreciated, it becomes useful 
to apply the developed tools of analysis of 2D spectra for 
more obscured molecular systems, for instance the FCP 
complex, and try to obtain more information about their 
structure / composition. 

Recently, a few multidimensional spectroscopy studies 
of FCP reported valuable insight about various unknown 
characteristics of the complex. A novel method of two- 
color two-dimensional spectroscopy allowed to map the 
the carotenoid-chlorophyll energy transfer pathways and 
resolved spectral heterogeneity of carotenoids in FCP 25 . 
It was found, that fucoxanthins in FCP are spectrally 
distinguishable and yet transfer energy to chlorophylls 
very efficiently. In our recent study^, dynamics of two- 
dimensional spectra of Chi a and Chi c Q y absorption 
bands at room temperature were analyzed and we re¬ 
ported ultrafast energy transfer between Chi c and Chi a, 
with signatures of excitonic coupling. However, limited 
temporal and spectral resolution in these particular ex¬ 
periments did not provide the possibilities to fully resolve 
coherent dynamics or a more fine structure of the FCP 
spectrum. 

Here we present results of the 2D ES of FCP obtained 
with higher spectral and temporal resolution than in our 
previous work 12 . To better discriminate the weak spec¬ 
tral features of FCP, experiments were carried out at 
77 K instead of the room temperature. The 2D spec¬ 
tra were analyzed with both femtosecond and picosecond 
time steps, thus, providing us the possibility to separate 
coherent and incoherent processes. In this work, we show 
that there are at least two Chi c molecules in FCP, which 
are clearly spectrally resolvable and participate in energy 
transfer on remarkably different timescales. We discuss 
that this might be related to certain spatial arrangement 
of Chi a and c molecules in the complex. Finally, the 
analysis of rich coherent dynamics during the first pi¬ 
cosecond reveals signatures of Chi a and fucoxanthin in¬ 
tramolecular vibrations, the former of which might play 
a role in enhancing ultrafast Chi c-to-a energy transfer. 


II. MATERIALS AND METHODS 

Fucoxanthin-protein complex was extracted and puri¬ 
fied from diatom Cyclotella meneghiniana. The proce¬ 
dure of purification and extraction is given in detail in 
Ref. 26 . The sample for the measurements was diluted 
with the glycero-buffer-detergent solution (10 mM Mes 
pH 6.5 buffer, 2 mM KC1, 0.03% (w/v) n-Dodecyl (3- 
D-maltoside and 60% (v/v) glycerol) to achieve optical 
density of about 0.3 at Chi a Q y band in 0.5 mm cell. 

The sample was placed in a 0.5 mm quartz cell and 
gradually cooled to 77 K in a liquid N 2 bath cryostat 
(Oxford Instruments). 

The setup for coherent two-dimensional electronic 
spectroscopy was based on the inherently phase- 
stabilized design by Brixner et alP3. In more detail, 
10 fs, 100 nm full width at half maximum (FWHM) 
light pulses (3 nJ energy per pulse at the sample) cen¬ 
tered at 615 nm were generated by home built non- 
collinear optical parametric amplifier were pumped by 
the Yb:KGW amplified laser system (Pharos, Light Con¬ 
version Ltd.) with repetition rate reduced to 7 kHz 
by means of pulse picker. Diffractive-optics-based non- 
collinear four-wave mixing setup with phase-matched box 
geometry was used. The third-order signal was obtained 
by spectral interferometry with the help of a heterodyn¬ 
ing local-oscillator pulse. Sensitivity and noise resistance 
was significantly enhanced by means of double modu¬ 
lation lock-in detection 28 ! Coherence time delay was 
scanned within -140-300 fs time range with 2 fs time step 
controlled by movable fused-silica wedges in beams k\ 
and & 2 - 2D spectra were obtained by processing the ac¬ 
cumulated sequence by the Fourier transform method 29 . 
The absolute phase of the 2D spectrum was retrieved 
by matching its projection with the transient absorption 
spectrum 30 recorded with the same setup with only two 
beams (for pump and probe) in use. Real (absorptive) 
parts of the total (rephasing+nonrephasing) spectra were 
used in further analysis. 

III. RESULTS 

The absorption spectrum of FCP measured at 77 K 
is shown in Fig. [l] by the black solid line. The spec¬ 
trum contains a well expressed band at 14895 cm -1 and a 
few weaker features around 15700 cm -1 and 17000 cm -1 . 
The dominating band at 14895 cm -1 is known to repre¬ 
sent the Q y transition of the core pigment Chi a. The 
feature at 15680 cm -1 was assigned to the Q y transition 
of Chi c, while the origin of bands at around 16200 cm -1 
and 17000 cm -1 were discussed to be reflecting the Q x 
transitions of Chi a and Chi c, respectivel}^. 

The spectrum of the laser pulse used in the 2D ES 
experiment is shown by a gray shaded area in Fig. [l] 
The spectrum was centered at around 16240 cm -1 with 
the FWHM of ^ 2600 cm -1 . The laser spectrum was 
intentionally shifted to the blue to enhance the signals 
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Figure 1. The absorption spectrum of FCP at 77 K (solid 
black line), square root of averaged (delay times from £2 — 
30 fs to 800 fs) diagonal cut of the 2D spectrum (dashed black 
line), excitation laser spectrum of 2D ES experiment (filled 
gray area). Thicker vertical green and blue sticks indicate 
the energies of the 0-0 Q y transitions of Chi a and Chi c, 
respectively. Vibrational states of Chi a, proposed by the hole 
burning measurements 3 ^ are shown by thinner green vertical 
sticks. Their heights indicate the Franck-Condon factors of 
the corresponding vibrations (values on the right axis). 


form the weak absorption bands of FCP: while it only 
partly covered the main absorption band at 14895 cm -1 , 
the other weak features were fully covered by the laser 
spectrum. 

Two different 2D ES experiments were carried out. In 
the first one, the 2D spectra at delay times £2 up to 800 fs 
were taken at every 10 fs. Such a high temporal resolu¬ 
tion is very useful for the analysis of the ultrafast ex¬ 
citation transfer dynamics as well as for studies of the 
quantum coherences and allowed us to resolve coherent 
beatings up to 1668 cm -1 in the frequency domain. The 
second set of data consisted of the 2D spectra measured 
at delay times £2 up to 1 ns with different time steps: 
from 50 fs for the initial evolution to 100 ps for delay 
times longer than 100 ps. These measurements were used 
for the analysis of the incoherent long-time population 
dynamics in FCP. 

A few 2D spectra at different delay times (£2 = 30 fs, 
330 fs, 5 ps and 30 ps) are depicted in Fig. [2K. In the 
spectra, the main band corresponding to the Chi a tran¬ 
sition at CJ3 = = 14926 cm -1 dominates. The fact 

that the band is not observed exactly at 14895 cm -1 as 
in the absorption spectrum, but rather blue-shifted by 
31cm -1 is caused by the excitation laser spectrum be¬ 
ing slightly off-resonant with respect to the main FCP 
absorption band (see Fig. TJ. The negative excited state 
absorption is well resolved above the diagonal in the 2D 
spectrum in the vicinity of the main positive absorption 
peak. The excited state absorption features in the 2DES 
could arise both due to excitonically coupled pigments 
(chlorophylls) and the corresponding double-excitation 
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as well as from the double-excitations of 
chlorophylls themselves 35 . 

In the 2D spectra shown in Fig. [2j thinner contour 
lines, drawn at every 2% of the maximum highlight 
the fine structure of the diagonal and off-diagonal peaks 
(cross-peaks) away from the main Chi a absorption band. 
The most pronounced features correspond to the cross¬ 
peaks lying around the same emission frequency CJ3 as 
the Chi a absorption band. These peaks might either 
signify the excitonic coupling between Chls a and c 12 { or 
they might be caused by the Chi a vibrational modes. 

In the region of the Q y of the Chi c transition fre¬ 
quency, i.e. at around CJ3 = u)\ ~ 15700 cm -1 , two di¬ 
agonal peaks are clearly resolved at excitation/emission 
frequencies 15610 cm -1 and 15810 cm -1 . The positions 
of these peaks are indicated by two vertical bars in the 
absorption spectrum shown in Fig. [I] coinciding with the 
maxima of the diagonal cut of the averaged 2D spec¬ 
trum (black dashed line in Fig. [l]). The zoomed-in re¬ 
gion of the 2D spectrum from CJ3 = uq = 15300 cm -1 
to CJ3 = uj\ = 16100 cm -1 (shown by a red square in 
Fig. [2]A.) is depicted in Fig. [2)3 at delay time £ 2 = 30 fs. 
No cross-peak associated with these two transitions is vis¬ 
ible at any delay time disproving the presence of excitonic 
coupling between states at 15610 cm -1 and 15810 cm -1 . 

The 2D spectra contain rich dynamics both at short 
(femtosecond) and long (picosecond) delay times, as it 
was extracted from the separate measurements described 
above. Dynamics of a few selected peaks in the spec¬ 
tra, occurring during first 30 ps are shown in Fig. [2p. 
The main diagonal peak (indicated by “1-1”) decays 
with timescales of 220 d= 100 fs and 5.4 db 0.6 ps. The 
cross-peaks (indicated by “2-1”, “3-1” and “4-1”) de¬ 
cay with the same longer rate of 5.4 ps as the “1-1” 
peak; the initial evolution is dominated by the high am¬ 
plitude coherences (Fig. [2p) and, thus, shorter lifetimes 
could not be unambiguously evaluated. Diagonal peaks 
at 15610 cm -1 and 15810 cm -1 show completely differ¬ 
ent long-time dynamics. The upper peak decays much 
slower than the lower one, as it can be seen by the 
spectral traces in Fig. [2p. We could fit the dynamics of 
both peaks as single-exponential decays with timescales 
of Ti56io = 320 d= 70 fs and rissio = 3.9 =b 1.6 ps, re¬ 
spectively. For longer delay times, additional exponential 
decays are present. 

The short-time (femtosecond) evolution of the same 
points in the spectrum are shown in Fig. HP demon¬ 
strating the oscillating behavior of the the diagonal and 
off-diagonal peaks corresponding to the Chi a. Coher¬ 
ence beatings can be analyzed by performing the Fourier 
transform of the time traces of the peaks. A very con¬ 
venient way to represent the overall complexity of os¬ 
cillations in the 2D spectrum is through the Fourier 
maps 36 41 . They are constructed by performing the 
Fourier transform of the each data point in the 2D spec¬ 
trum over the delay time £2 after subtraction of slow ex¬ 
ponential decay contributions. Then, the Fourier am¬ 
plitude can be plotted as a two-dimensional map against 
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Figure 2. (A) 2D absorptive spectra at different delay times. All the spectra are normalized to the maximum of the spectrum 
at £2 = 30 fs and drawn using the same scale. Thicker contour lines are drawn at every 10% of the maximum of the scale; 
thinner lines are drawn every 2%, up to ±20%. (B) The two-peak structure in the 2D spectrum of FCP. Zoomed-in region, as 
shown in Fig. [2j\ by the red square. Values in the spectrum are multiplied by 20 and the spectrum is drawn using the same 
scale as in (A). (C) Picosecond and (D) femtosecond time traces of a few selected peaks indicated by arrows in (A) and (B). 
Values obtained by averaging around the selected points within the circle of 100 cm -1 radius are presented. Multi-exponential 
fits of dynamics of a few peaks are shown by the blue dashed lines (see the text for the fitting parameters). 


excitation and detection frequencies cand CJ3 for a fixed 
Fourier frequency Each Fourier map then can be inte¬ 
grated over lji and UJ 3 and the integral value represented 
as function of coherence frequency <^2- Such a depen¬ 
dency would show the averaged distribution of various 
oscillation frequencies. 

We chose to represent the amplitude of observed oscil¬ 
lations by the value of the Fourier map norm, which is 
equal to the square root of the sum of squares of the data 
points of Fourier maps at frequency u) 2 - From the oscil¬ 
lation amplitude dependency shown in Fig. [3]A, 5 sharp 
peaks at = 255 cm -1 , 345 cm -1 , 745 cm -1 1155 cm -1 
and 1515 cm -1 can be readily pointed out. More rigor¬ 
ous analysis of the second derivative of the oscillation 
amplitude function reveals 24 more frequencies. 


IV. DISCUSSION 

The 2D ES spectra obtained with femtosecond and pi¬ 
cosecond time resolutions give a wealth of information re¬ 
garding both evolution of population and coherent states 
of the molecular system under consideration. Special at¬ 
tention should be taken into consideration when analyz¬ 
ing the processes occurring on different timescales. Thus, 
before turning to the analysis of the femtosecond evolu¬ 
tion of the measured spectra, let us discuss a few mech¬ 
anisms that contribute to the spectral lineshape forma¬ 
tion: pulse overlap, spectral diffusion, dynamic Stokes 
shift and coherence beatings. They can be important 
factors that might lead to ambiguities and errors when 
interpreting the time-dependencies of the 2D spectrum. 

Pulse overlap effects dominate during the initial de- 
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Figure 3. (A) Amplitude of oscillations (Fourier map norm) as 
a function of coherence frequency Dependency extracted 
from Fourier transforms of zero-padded and not-padded spec¬ 
tral time evolutions are shown by symbols and lines, respec¬ 
tively. Determined frequencies of oscillations are shown by 
the vertical black and gray lines. (B) Franck-Condon factors 
of vibrational intramolecular modes of Chi a from the hole 
burning experiment 31 (green sticks), proposed frequencies of 
electronic coherences (blue dashed lines) and vibrational Ra¬ 
man frequencies of fucoxanthin (orange dashed lines). 


lay time evolution, comparable to the length of the laser 
pulses. That is, besides the double-coherence path¬ 
ways, additional rephasing contributions appear in the 
phase matching direction when the k\ and pulses do 
overlap 42 . Therefore, the initial evolution of the 2D spec¬ 
trum (in our case first 30 fs) is dominated by the sharp 
decay of these additional pulse-overlap contributions and 
should be considered very carefully as they could be con¬ 
fused with some ultra-fast relaxation mechanisms. 

Spectral diffusion can also make an evident influence 
on the 2D spectrumPSl. It accounts for the loss of correla¬ 
tion of the narrow-bandwidth excitation and spectrum is 
manifested by the increase of the anti-diagonal width of 
the peaks in the absorptive 2D spectrum with the delay 
time. The spectral diffusion affects both the rephasing 
and nonrephasing spectral contributions. For a simple 
two-level system the rate of the spectral diffusion can 
be related to the bath relaxation time and results in ex¬ 
ponential decrease of the height of the main absorption 
peak in the 2D spectrunP^. As such an exponential evo¬ 
lution overlaps with the other temporal lineshape forma¬ 
tion mechanisms, it becomes troublesome not to over¬ 
interpret the spectral diffusion as a relaxation-type sig¬ 
nal. Comparing the evolution of the peak height against 
its anti-diagonal width provides a possible way to evalu¬ 
ate the contribution of the spectral diffusion. However, 
sometimies this is not possible because of overlaps of 
closely-positioned coupled states and cross-peaks. The 
other method is to check the time evolution of the inte¬ 
gral over some area of the spectrum in the vicinity of the 
peak under consideration in the absorptive spectrum. If 
the radius of the integrated area is larger than the homo¬ 
geneous linewidth, the spectral diffusion is averaged out 
from such signal. 

System states during the delay time £2 experience re¬ 
laxation processes through the vibronic pathways. Such 


a vibrational cooling is often referred to as the dynamic 
Stokes shift and expressing itself as the shifting of peaks 
down in the energy scale along the emission axis with £2 
in the 2D spectrum. As the ground state bleaching and 
excited state emission contributions can, in principle, ex¬ 
perience the dynamic Stokes shift on different timescales, 
it can cause some overlaps of these positive and negative 
contributions, giving a complex time-dependence of both 
diagonal peaks and cross-peaks. Moreover, the timescale 
of the dynamic Stokes shift can be around a few hundreds 
of femtoseconds - therefore, without a careful analysis 
such evolution might once again be wrongly assigned to 
some ultrafast relaxation pathway. 

Finally, the initial evolution of the 2D spectrum is se¬ 
riously affected by the coherence beatings. Lifetimes of 
the electronic coherences can be as short as tens of fem¬ 
toseconds and vibrational/vibronic coherences can result 
in oscillations for as long as a few picoseconds. Although 
static disorder tends to dephase electronic coherences 
very quickly, their contribution to the initial dynamics 
can still be substantial. For example, in spectra recorded 
with a low £2 resolution, the coherent signal adding on 
top of the population evolution can be as significant as 
the population signal itself. Thus, coherence beatings 
can seriously impede interpretations of the population 
dynamics if the 2D spectrum is of low resolution in delay 
time £2- 


A. Identification of coherences 

In the 2D ES spectra of FCP we would expect co¬ 
herences of various origins and amplitudes to appear: 
ground state vibrational coherences of Chi a, Chi c and 
Fx, excited state vibronic coherences of the same fre¬ 
quencies, mixed coherences of shifted or combinational 
frequencie d 16 * 45 * and electronic coherences due to exci- 
tonic coupling of chlorophylls. Such distinction between 
the coherences of different origin is conveniently observed 
in the Fourier amplitude maps. In the rephasing spec¬ 
trum, ground state coherences are manifested by strong 
beatings below and on the diagonal 38 while the corre- 
spondingvibronic coherences - symmetrically above the 
diagonaP®^. Mixed coherences, that represent super¬ 
positions of quantum mechanically mixed states, are ob¬ 
served as complex and as ymme tric patterns of oscillating 
peaks in the Fourier mapd^EU. Electronic coherences in 
the rephasing spectrum manifest themselves exclusively 
as oscillatory cross-peaks and no oscillations on the di¬ 
agonal and vice versa for the nonrephasing spectrum 38 . 
However, oscillations in the nonrephasing contribution 
of the spectrum are more sensitive to static disorder of 
pigment energies and, thus, in the systems with large 
inhomogeneous broadening coherences from the rephas¬ 
ing contributions will dominate over those from the non¬ 
rephasing contribution. 

Let us discuss about the origin of coherences in the 
spectra of FCP by comparing the obtained dependencies 
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of the intensity of oscillations (Fig. [3]4) with known fre¬ 
quencies corresponding to the intramolecular vibrational 
modes of Chi a . We refer to the hole burning exper¬ 
iments on single chlorophylls carried out in 1.6 K by 
Gillie et alP^ I (green sticks in Fig. [3^) and a recent study 
of Du et alP^, in which the frequencies were obtained 
from the time and frequency resolved pump-probe spec¬ 
tra at 297 K by the Fourier analysis, similar to the one 
employed here. 

Frequencies, reported in these studies as well as ob¬ 
tained by our experiment, are summarized in Table [l| 
The comparison gives almost perfect agreement not only 
for the obtained strongest coherences at around UJ 2 = 
255 cm -1 , 345 cm -1 , and 745 cm -1 (deviations smaller 
than 11cm -1 ) but also for the weaker ones with the 
largest deviation of 20 cm -1 for the UJ 2 = 1305 cm -1 vi¬ 
brational mode. The Fourier amplitude maps at co >2 = 
255 cm -1 , 345 cm -1 , and 745 cm -1 are shown in Fig. [ 4 ] 
All three of them clearly indicate the presence of vibra¬ 
tional and vibronic coherences as the patterns of the oscil¬ 
latory features can be easily decomposed to 4 stimulated 
emission (SE) and 4 ground state bleaching (GSB) contri¬ 
butions of a simple two-level system with one vibrational 
state in the excited and ground states (see the scheme 
and the corresponding excitation pathways in Fig. 
and C). From the scheme it follows that weak features 
above the diagonal, separated by exactly UO 2 value, indi¬ 
cate vibronic excited state coherences (SEi pathway). 

The map at UJ 2 = 745 cm -1 lacks amplitude of oscil¬ 
lations on the diagonal compared to the UJ 2 = 345 cm -1 
or UJ 2 = 255 cm -1 maps due to the narrow bandwidth of 
the laser spectrum, suppressing some specific coherences. 
That is, since Chi a 0-0 transition is on the low-energy 
shoulder of the laser excitation spectrum, all excitation 
pathways containing lower energy transitions (indicated 
by the dashed arrows in Fig. |5p) will be substantially 
suppressed. These particular pathways in Fig. HP are 
shown in a shaded background. 

For oscillations of higher frequencies, the mentioned 
contributions will be completely excluded. Therefore, 
the dependency of amplitude of oscillations in Fig. [3]A. 
should be considered carefully, as the peak heights there 
do not perfectly represent the amplitudes of all summed 
pathways, but are rather suppressed in middle and high 
frequencies. 

The Fourier amplitude maps at u )2 = 1155 cm -1 
and 1515 cm -1 are shown in Fig. [ 4 ] They both con¬ 
tain broad features in the region of excitation frequency 
uji = 16500 — 17000 cm -1 and emission frequency CJ 3 ~ 
uj\ — 1155 cm -1 and CJ 3 ~ uq — 1515 cm -1 , respectively. 
With our laser spectrum we could in principle excite 
the S 2 states of some of the “red” fucoxanthins and, 
possibly, get signatures of carotenoid ground state co¬ 
herences (from the GSBf* pathway shown in Fig. [HJl)). 
These contributions would appear as oscillations at the 
emission frequency, shifted down from the diagonal by 
the frequency of the vibrational mode. Therefore, we 
assign these features observed at UJ 2 = 1155 cm -1 and 
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Table I. Summary of coherent oscillations (cm -1 ) in the mea¬ 
sured FCP 2D spectrum (first column). Vibrational frequen¬ 
cies of chlorophyll a obtained from the hole burning 3 ^ (HB, 
second column) and pump-probe^ (PP, third row) experi¬ 
ments. Vibrational frequencies of fucoxanthin from the reso¬ 
nant Raman (RR) study^ are shown in the fourth row. Num¬ 
bers in bold denote the five most intensive oscillations in 
2D and pump-probe experiments, two most intensive Raman 
modes of fucoxanthin and frequencies from the hole-burning 
experiment, having Franck-Condon factors higher than 0.02. 


UJ 2 = 1515 cm 1 to the ground state vibrational modes 
of Fx at 1160 cm -1 and 1536 cm -1 , respectively. 

Apart from the Fx vibrational coherence, the Fourier 
amplitude map at UJ 2 = 1515 cm -1 shows oscillations at 
CJ 3 = 14935 cm -1 and uq = 16550 cm -1 and a corre¬ 
sponding diagonal-symmetric feature. They can be as¬ 
signed to the Chi a vibrational and vibronic coherences, 
respectively. Thus, the full pattern of features in the map 
at CJ 2 p= 1515 cm -1 can be explained by the scheme shown 
in Fig. [5jl: the two diagonal-symmetric narrow features 
come from the excited state vibronic and ground state 
vibrational coherences of Chi a and the broad feature 
below the diagonal is due to vibrational ground state co¬ 
herences of Fx. The other weak features at the emission 
frequency CJ 3 ~ 16000 cm -1 cannot be explained by any 
sort of excitation pathways; we assign them to experi¬ 
mental artifacts coming from the nonrephasing contribu¬ 
tion (not shown). 

The Fourier map at U 02 = 1155 cm -1 clearly lacks sim¬ 
ilar features stemming from Chi a vibrational coherences 
and there is no vibrational mode of similar frequency sug¬ 
gested by the hole-burning or pump-probe experiments 
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Figure 4. Fourier amplitude maps at coherence frequencies UJ 2 — 255 cm -1 , 345 cm -1 , 745 cm -1 , 1155 cm -1 and 1515 cm -1 . 
Fourier amplitude is shown by a color scale, normalized to the maximum of the map at U 2 — 345 cm -1 . Black contour lines 
show the 2D spectrum at £2 = 30 fs. Diagonal and parallel to the diagonal dashed lines are shifted by the corresponding value 
of oj 2 . Notice that excitation axex of maps at 0 J 2 = 1155 cm -1 and 1515 cm -1 are extended compared to the other maps. 


0 J 2 — 345 cm 1 B CJ 2 — 1515 cm 1 






Figure 5. Patterns of features in the Fourier amplitude maps 
at UJ 2 — 345 cm -1 (A) and UJ 2 — 1515 cm -1 (B). Excitation 
pathways, leading to coherences in these maps are shown in 

(C) and (D), respectively. Horizontal thick lines in (C) and 

(D) show ground and excited state energies of Chi a and Fx. 
Thin horizontal lines indicate the corresponding vibrational 
states. Vertical straight arrows represent interactions; verti¬ 
cal wavy lines indicate emission. Transitions, suppressed by 
the laser spectrum are shown by dashed lines and the path¬ 
ways, becoming negligible for the higher frequency vibrational 
modes, are shown in the shaded background. 


considered here 3135 . However, vibrational Chi a mode 
of 1145 cm -1 has been reported in the fluorescence nar¬ 
rowing and resonant Raman study by Telfer et al 49 ( It 
is possible that we were not able to resolve signatures of 
this particular vibrational coherence because of the sup¬ 
pression of high frequency coherences by laser spectrum 
as discussed above. Alternatively, there could be addi¬ 


tional mechanisms, responsible for some vibrational Chi a 
frequencies not being detected in FCP, e.g. vibronic cou¬ 
pling between Chi a in FCP, adding more complexity to 
spectra, similar to the enhancement of the ground state 
coherence proposed by Chenu et aP^, or the anharmonic- 
ity of the electron-vibrational coupling recently proposed 
for the LH2 antenad^. Obviously, the situation calls for 
more delaited study of vibrational/vibronic coherences of 
chlorophylls in coupled pigment-protein complexes. 

After performing analysis of the Fourier maps of the 
other coherences, we were not able to extract any os¬ 
cillatory features that could be assigned to the vibra¬ 
tional/vibronic coherences of Chi c. This is probably due 
to too weak signal. Indeed, the Chi c absorption signal 
is at least 5 times weaker than the main absorption band 
of Chi a. Thus, the signal of coherent oscillations might 
become too weak to be distinguished from the exper¬ 
imental noise. A close inspection of oscillations around 
the proposed Q y transition band of Chi c 15600 cm -1 ) 
revealed only very weak signatures of vibrational coher¬ 
ences of Chi a and Fx. The former might be coming from 
the nonrephasing vibrational contribution and the latter 
- from the overlap of a wide peak of the fucoxanthin 
ground state vibrational coherence. 

Obviously, the strongest high-frequency electronic co¬ 
herence in the spectra of FCP would be the coherence of 
Q y transitions of excitonically coupled Chi c and Chi a. 
From the position of the corresponding peaks in the ab¬ 
sorption spectrum, the frequency of such a coherence 
should be around 715 cm -1 or 915 cm -1 . By comparing 
the peak heights and assuming that the Franck-Condon 
factors of vibrational transitions are as small as 0.05, it 
could easily be deduced that the amplitude of electronic 
coherence should be higher than the vibrational coher¬ 
ence of Chi c. But the dephasing rate of the electronic 
coherences is known to be much shorter than the vibra¬ 
tional oned^HI. This would prevent the electronic coher¬ 
ence to be captured by the Fourier maps, except of a few 
cases. One of them is the inhomogeneous distribution in¬ 
duced picosecond electronic coherence lifetime, as it was 
recently shown by Dong et alPl. In the inhomogeneously 
broadened spectrum the sub-ensemble of molecules, hav- 







































































































































ing similar energies provides long-lived coherent signal 
in the two-color photon echo experiment. Such mecha¬ 
nism could also be significant in the 2D electronic spec¬ 
troscopy, however, it requires the absorption spectra of 
Chi c and Chi a to overlap significantly. 

Another mechanism allowing for long-lived electronic 
coherences is substantial electronic-vibrational mixing 
in the excited state, resulting in a vibrational coher¬ 
ence lifetime borrowing mechanism 47 ! For this effect 
to be noticeable, the excitonic energy gap between the 
excitonic states should closely match some strong dis¬ 
crete vibrational mode. In our case, the proposed ex¬ 
citonic energy gap is close to a very strong vibrational 
modes of 745 cm -1 or 915 cm -1 and in the Fourier am¬ 
plitude maps the ground state vibrational contribution 
dominates. Therefore, if there are observable signatures 
of electronic coherence beatings, they are hindered by 
strong signals of vibrational coherences. 


B. Ultrafast Chi c-Chl a transfer 

In our previous study of the 2 D spectra of the FCP 
complex at room temperature, the ultrafast energy trans¬ 
fer (occurring on timescale of 60 fs) from excitonically 
coupled Chi c to Chi a was proposed 12 . It was obtained 
by tracing the time dynamics of the spectral feature at 
CJ 3 = = 15680 cm -1 , assigned to the Q y transition 

of Chi c. Here we report that there are two peaks at 
15610 cm -1 and 15810 cm -1 instead of one. This doublet 
was not resolved in the experiment at the room tem¬ 
perature. Moreover, we found significant differences in 
their decay dynamics. Considering times up to 30 ps, 
the lower peak decays with the timescale of 320 fem¬ 
toseconds, while the the shortest temporal dynamics of 
the upper peak is 3.9 ps. The timescales were obtained 
from temporal traces, taken by averaging over the area of 
a circle of 100 cm -1 radius in spectra and omitting the 
first 30 fs, thus limiting amount spectral diffusion and 
pulse overlap contributions appearing in peak dynamics. 
We can thus conclude that the obtained two different 
Chi c states directly transfer energy to Chi a on different 
timescales. Notice that energy transfer via intermediate 
states (presumably from fucoxanthins) is not consistent 
with fluorescence-excitation spectrum of FCP^, which 
implies 100% efficiency of Chi c to Chi a energy transfer, 
and the fact that significant portions of excitation from 
Fx Si/ICT (internal charge transfer) state, which could 
be suggested to be a bridge between Chi c and Chi a, are 
lost due to the decay to the ground state 25 . 

Recent studies suggest that there are two Chi c 
molecules in the comple^. Also, some of them can be of 
different types (i.e. Chi ci, C 2 or C 3 ) that have slightly dif¬ 
ferent absorption spectra 53 . Even if the discussed spec¬ 
tral features are caused by the chlorophylls of different 
types (thus, causing signals in different frequencies), the 
temporal dynamics of the corresponding peaks indicate 
that their interaction with the rest of the system or a 


surrounding environment is different. The simplest pic¬ 
ture is that one Chi c is situated in the periphery of 
the FCP complex and is weakly coupled to the rest of 
Chi a, thus, the excitation transfer to the core pigments 
is slow and we observe a picosecond decay of the peak at 
15810 cm -1 . At the same time, the protein environment 
causes it to have higher energy. Another Chi c, absorb¬ 
ing at 15610 cm -1 , must be coupled to Chi a stronger for 
faster excitation transfer. Since the cross-peak between 
these two transitions as well as the relevant excited state 
absorption contribution is absent, chlorophylls c in FCP 
are not coupled. This implies that they probably are 
spatially well separated. 

However, the ultrafast energy transfer from the lower 
energy Chi c state still remains an issue to be explained. 
The strength of coupling between Chi c and Chi a is not 
known, although it may be estimated from LHCII, which 
shares considerable sequence homolog}®. Couplings be¬ 
tween chlorophylls in the LHCII complex were calculated 
to be less than 100 cm - ®. Considering this and the 
fact that the transition dipole moment of the Q y band 
is smaller for Chi c than for Chi a, one can be certain 
that coupling between Chi c and Chi a in FCP should 
be even weaker. This is also supported by the absence 
of clear excitonic cross-peaks between Chi a and Chi c 
states in our 2D spectra. 

The Forster energy transfer could apply in this range 
of parameters. However, it leads to slower picosecond 
transfer rates. The 320 fs transfer timescale observed in 
2 D ES can be understood by including two ingredients 
in the consideration. First, the lower energy Chi c can 
be coupled to several Chi a molecules, thus, it transfers 
energy to a few Chi a molecules at once. The rates ef¬ 
fectively sum-up and lead to the fast decay of the Chi c 
population. Second, vibrations-assisted energy transfer 
enhancement might be significant for the Chi c-Chl a 
energy transfer. The energy gap of 715 cm -1 is close to 
the strongest vibrational Chi a mode of 745 cm -1 , thus 
allowing for electronic-vibrational resonanc^^P Re¬ 
cently, specific vibrational modes were suggested to speed 
up the initial charge tra nsfe r step in the photosystem II 
(PS II) reaction centei^^. Finally, this does not con¬ 
tradict our previous suggestion that Chi c is coherently 
coupled to Chi cP®, however, the coherent mechanism 
may be far more complicated and calls for future studies. 

We propose a schematic spatial arrangement of Chi a 
and Chi c molecules (shown in Fig. [ 6 |. This functional 
scheme is comparable to the CP29 photosynthetic an¬ 
tenna of photosystem l l 58 * 59 i The structure of the CP29 
complex is similar to the light harvesting complex LHCII, 
which shares the sequence homology with FCP. In CP29 
there are eight Chi a and four Chi b (the absorption maxi¬ 
mum of Chi b is around 15674 cm -1 , similarly as to Chi c) 
and one more chlorophyll binding site can be occupied by 
either type of chlorophyll, although previously it was 
assumed that there are 8 Chi a and 2 Chi b molecule^. 
However, in both models one chlorophyll b molecule is 
well spatially separated from the other which is strongly 
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Figure 6. Proposed functional scheme of spatial arrangement 
of chlorophyll c and a molecules in the FCP complex. Cou¬ 
plings between pigments are shown by solid lines. Dashed 
lines indicate weaker coupling between spatially separated 
Chi c absorbing at 15810 cm -1 . Chi a molecules shown here 
do not represent their actual number in FCP. 


coupled to the rest of chlorophylls a of the complex. Con¬ 
sequentially, two absorption bands at 15650 cm -1 and 
15385 cm -1 are observed in the absorption and in the 
pump-probe spectrum have completely different tempo¬ 
ral dynamics: the energy transfer from Chi b at higher 
energy is transferred to Chi a with the rate of 350 fs, 
while transfer from the lower energy Chi b band occurs in 
2.2 ps 58 . This implies that having two accessory chloro¬ 
phylls that are in different spatial situations might be 
a property shared by a number of light-harvesting com¬ 
plexes, thus probably it has a functional significance. 

In a recent study of FCP based on resonance Raman 
techniqu^, signatures of two distinct Chi c molecules 
were observed. In that study a possible pigment organi¬ 
zation scheme in FCP was discussed, based on the pig¬ 
ment binding sites of the LHCII from plants 10 . It was 
suggested that one of the Chi c should be bound at either 
the site of a614 (which is preserved in FCP) or a613 (not 
preserved in FCP), while the other Chi c should be bound 
at the site of b609 (preserved) or a604 (not preserved). 
For the faster transferring Chi c the site of b609 would 
be a logical choice, as it is in the vicinity of several pig¬ 
ments. Meanwhile, the higher energy (15810 cm -1 ) Chi c 
could indeed be situated at the site of a614, as that site 
is relatively isolated and would be in accord with slower 
energy transfer from this pigment. Since this pigment 
is situated at the periphery of the complex, it might be 
responsible for the inter-complex Chi c to Chi a transfer 
in the FCP-PSII supercomplexed^. 


V. CONCLUSIONS 

In this work we studied low-temperature and high- 
resolution experimental two-dimensional spectra of FCP 
complex. Short time (femtosecond) dynamics of the spec¬ 


tra show complex oscillation patterns. They are dom¬ 
inated by vibrational and vibronic coherences of Chi a 
and fucoxanthin; electronic coherences could not be iden¬ 
tified due to their short lifetime and being hindered under 
high-amplitude oscillations originating from Chi a. Most 
notably, we have resolved two distinct states, that can be 
assigned to chlorophylls c. These states were not found 
in the previous studies of FCP and their energies were 
not known. The two Chi c molecules transfer energy to 
Chi a with rates that differ by about an order of magni¬ 
tude. We suggest that the lower energy Chi c molecule 
is situated in the core of the complex and is coupled to 
an aggregated pool of Chi a molecules and the Chi c-to- 
a energy transfer might be vibronically-enhanced. The 
higher energy Chi c should be in the periphery of the 
complex to transfer energy on the picosecond timescale. 
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